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ABSTRACT 

We demonstrate in this paper that the velocity widths of the neutral gas in Damped Lya (DLA) systems are 
inconsistent with these systems originating in gas disks of galaxies similar to those seen in the local universe. 
We examine the gas kinematics of local galaxies using the high quality H I 21 -cm data from the H I Nearby 
' Galaxies Survey (THINGS) and make a comparison with the velocity profiles measured in the low-ionization 

' metal lines observed in DLAs at high redshifts. The median velocity width of z = H I gas above the DLA 

, column density limit of A'hi = 2 x 10^" cm^^ is approximately 30 km s^^, whereas the typical value in DLAs 

C\| ■ is a factor of two higher. We argue that the gas kinematics at higher redshifts are increasingly influenced by gas 

^ that is not participating in ordered rotation in cold disks, but is more likely associated with tidal gas related to 

O , galaxy interactions or processes such as superwinds and outflows. An analysis of the H I in the local interacting 

■ star-burst galaxy M82 shows that the velocity widths in this galaxy are indeed similar to what is seen in DLAs. 

, Subject headings: galaxies: ISM — ISM: evolution — ISM: kinematics and dynamics — radio lines: ISM — 

galaxies: statistics — quasars: absorption lines 
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1. EVITRODUCTION 

Until the next generation of large radio telescopes system- 
atically detects the 21 -cm hyperfine line of neutral, atomic 
hydrogen (Hi) in emission at redshifts beyond z > 0.2, our 
knowledge of the H I at high redshift will remain solely de- 
pendent on the study of absorption-line spectra of high red- 
shift quasars. These spectra are dominated by Lya absorption 
lines originating in intervening H I clouds, and show a large 
range in Hi column density. The systems with the largest 
column densities, termed Damped Lya Absorbers (DLAs; 
IWolfe et al. 1986), are characterized by iVni > 2 x 10^" cm"^ 
and are known to hold most of the H I atoms in the uni verse, 
at all redshifts (cf. lPeroux et"ani2005l:IZwaan et al.ll2005 b). A 
popular explanation is that the DLAs are the prog enitors of 
gas disks in present-day galaxies ( Wolfe et"ani2005h . 

At z = it is well established that most of the H I atoms are 
locked up in the gas disks of _L,-type galaxies (IZwaan et al.l 
l2003h . These z = Hi disks show an frequency dis- 
tribution of Hi column densities [/(iVHi )1 that is virtually 
identical in shape to that of the DLAs (iRvan- Weber et al.l 
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2005h . Further^ 



12003; Zwaan et al. 2005bt iProchaska et"al 
more, at redshifts z < 1 there are approximately 20 identifi 
cations of galaxies gi ving rise to DL A absorption in spectra of 
background quasars ( Le Brun et al. 1997; Steidel et al. 1995 
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Lacy et ani2003|) . The properties of these galaxies (such as 
luminosities and impact parameters) are fully consistent with 
the idea that the DLAs arise in galax ies that are typical in 
the local universe (IZwaan et al.ll2005"bh . Also the metallicities 
seen in low redshift DLAs are fully consistent with what is 
expected when gas disks of z = galaxies are intercepted 
randomly. Therefore, DLAs at low redshifts (z < 1) proba- 
bly arise in gas disks of galaxies like those seen in the local 
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At higher redshifts (z > 1) the situation becomes more 
complicated. Direct detection of DLA host galaxies has 
prove n to be extremely difficult (see e.g.. lM0ller et al.ll2002l 
12004 and refer ences therein). This may not be surprising if at 
these epochs the DLAs also arise in gas disks of normal spi- 
ral galaxies: the cross-section selection ensures that a large 
fraction of the DLA galaxies is fainter than L,, and thus dif- 
ficult to detect (Fvnb o et al.l ll999; Zwa an et al.li2005bi) . The 
shape of the column density distribution function, /(A'hi), at 
these redshifts is still consistent with what is seen at z = 0. 
Although this is consistent with the galaxy disk model, it is 
inconclusive evidence for this model. It is conceivable that 
the shape of /(A^hi) is determined by rather simple physics, 
for example ionization at low column densities and H2 forma- 
tion at high densities. This might produce the same /(A'hi) 
at all redshifts, regardless of whether this gas is distributed in 
ordered disks or not. 

Ot her observations o f high-z DLAs are also inconclusive 
(see IWolfe et al.l I2005L for a comprehensive review). The 
inferred cosmic star formation rate density from DLAs is 
consis tent with th at measured from Lyman Break Galaxies 
(Wolf ed al.ll2003h . However, m etallicities of the DLAs are 
typically low (e.g., Pettini et al. 1999), as are the molecular 
fractions (e.g., L edoux et al. 2003) , and analyses of elemental 
abundances show s tar formation histories typical of those of 
quiescent galaxies (iDessauges-Zavadsky et al.ll2006h . 

Perhaps one of the most conclusive tests of whether DLAs 
at high redshift arise in gas disks or not, can be made by study- 
ing their kinematics. The first wor k on this topic was pre- 
sented by IProchaska & Wolfd (1 19971) . who compiled high res- 
olution, high signal-to-noise spectra of low-ionisation metal 
lines in DLAs, and compared them with the velocity pro- 
files derived from simulations using model ensembles of ar- 
tifical clouds with a range of velocities. They concluded that 
rapidly rotating thick disks are th e most likely ho s ts of D LAs. 
This conclusion was contested bv' Haehnelt et al.l(ll998b . who 
showed that irregular protogalactic clumps can reproduce the 
absorption line properties of DLAs equally well. A direct 
co mparison with 21 -cm d ata of local galaxies was first made 
by IProchaska et al.l (l2002l) . who found that low mass galax- 
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ies such as the Large Magellanic Cloud (LMC) display much 
lower velocity widths than normally seen in DLAs. 

A fair comparison with DLA kinematics can only be made 
by studying a representative sample of local galaxies for 
which high quality 21 -cm data cubes exist, and which have 
high spatial and frequency resolution and low H I column den- 
sity limits. In this paper we make use of the THINGS (The 
H I Neai-by Galaxy Survev: [Walter et al.l2008h data, which are 
ideal for this purpose. 

It is important to note that the Lya lines that define the 
DLA carry little information on the kinematics of the ab- 
sorbing systems themselves. The velocity width that would 
correspond to the total extent of the damping wings is much 
larger than the physical velocity width corresponding to the 
turbulence or large scale movements of the gas; any struc- 
ture in the velocity distribution of the Hi gas is completely 
washed out by the damping wings. Instead, workers in this 
field rely on low-ionization metal lines, such as Fe+, Si+, 
Ni"*", which are believed to be g ood tracers of the neutral gas 
(e.g.. Prochaska & Wolfdll997l) . In the same system each of 
these transitions gives consistent velocity widths, as long as 
the lines are unsaturated. In principle, many small clouds in 
the outskirts of galaxy halos could contribute to the total ve- 
locity profile measured in these low-ionization lines. How- 
ever, the parameter that we use to characterize the velocity 
width, AVgo (the velocity width that encompasses the central 
90% of the total optical depth seen in the line) ensures that 
these low column density clouds do not contribute to the total 
width. 

In section|2]we describe our strategy and data, in section[3] 
we present the results, and in section |4] we discuss these re- 
sults. Finally, we present the conclusions section|5] Through- 
out this paper we use a value for the Hubble constant of 
Hq = 75kms~^ Mpc^^, although the distances to some of 
the THINGS galaxies are derived using direct distance mea- 
surements. 

2. THE APPROACH 

For a quantitative comparison between the kinematics of 
DLAs and local galaxies, high quality 21 -cm maps are re- 
quired that are characterized by high signal to noise and high 
spatial and velocity resolution. A low noise level is required 
to sample all regions of the gas disks above the DLA column 
density limit of Nm — 2 x 10^°cm^^. A high spatial reso- 
lution is needed to avoid severe overestim ation of the veloc- 
ity pr ofiles due to spatial smoothing (see iRvan- Weber et al.l 
l2005h . Since for aperture synthesis interferometry the noise 
level (measured as the minimal detectable H I column density 
per beam) scales with the inverse of the beam area, these two 
requirements are very difficult to meet. The THINGS sample 
is currently the only sample with sufficient spatial resolution 
and depth, and contains a large enough number of galaxies 
such that the local universe galaxy population can be sampled 
representatively. 

2.1. THINGS 

We base our analysis on VLA observations of a sample of 
nearby galaxies drawn from THINGS: The HI Nearby Galaxy 
Survey (I Walter et al.ll2008h . THINGS provides high angu- 
lar and velocity resolution data of uniform quality for 34 
galaxies spanning a range of Hubble types, from dwarf ir- 
regular galaxies to massive spirals. The sample selection is 
mainly draw n from the Spitzer Inf rared Nearby Galaxy Sur- 
vey (SINGS. iKennicutt et ani2003b but excludes targets with 



distances > 15 Mpc, and elliptical galaxies. In addition, edge- 
on galaxies have been excluded. The observations an d data re- 
duction are discussed in detail in'W alter et alj (1200 8). For our 
analysis, we use the robust-weighted data cubes, which pro- 
duce typical beam sizes of 6" (or 200 pc for the mean distance 
of 6.5 Mpc to our sample galaxies). We choose the robust- 
weighted data cubes in order to achieve the highest possi- 
ble spatial resolution, which is required for making a com- 
parison to DLA velocity profiles, which are measured along 
very narrow sight-lines to background quasars. The veloc- 
ity resolution in all cases is equal or better than 5.2 kms^^ 
(which we will smooth further, see below). Typical column 
density sensitivities are 5 x 10^^ cm^^ in the integrated Hi 
maps/profiles, i.e., much below the DLA column density limit 
(2 X 10^° cm^^) discussed in this paper 

The different distances among the galaxies in the sample 
cause the beam to trace different physical scales within galax- 
ies. As discussed in section |231 beam smearing causes some 
increase in the measured velocity widths. Therefore, for the 
most distant galaxies the velocity width measurements are in- 
creased somewhat more than for the nearby ones. However, 
as we will see, this does not affect our main conclusions. 

2.2. Biases in the THINGS sample 

The THINGS sample is composed of galaxies from opti- 
cal c atalogues, wh ereas DLA systems are purely H l-selected. 
(see [Walter et al .''2008, for a discu ssion on the selection of 
THINGS targets). However, in Zwaan et all ( |2005b) it was ar- 
gued that this optical selection does not lead to any bi ases for 
this application. The main argument is the following: I ZwaanI 
(2000) has shown that the properties of Hl-selected galaxies 
are not different from those of optically selected galaxies if 
compared at the same luminosity. Therefore, if the weighting 
of the sample is properly taken into account (using the lumi- 
nosity function or H I mass function) no strong biases should 
be introduced by using optically selected galaxi es. This fact 
is also corroborated by the calculation of R yan- Weber et al.l 
12003), who showed that the redshift number density dN/dz 
calculated from H l-sel ected nearby galax ies is consistent with 
the values presented in lZwaan et al.l(l2005h) and in the present 
paper 

Another consideration is the distribution of inclinations i 
of the THINGS sample. Ideally, the galaxies should be ran- 
domly oriented, that is, the distribution of cos i should be con- 
sistent with a flat distribution. In practice, however, given 
the selection discussed in section 12.11 the sample has a mi- 
nor deficiency of edge-on galaxies. A Kolmogorov-Smirnov 
test shows that the null-hypothesis that the inclination distri- 
bution is consistent with being random can only be rejected 
at the 19% level. Because of projection effects, the most 
inclined galaxies have H I column densities and AVgo mea- 
surements at the high end of the total distribution of these 
values in the THINGS sample. On the other hand, the cross- 
section in highly inclined galaxies is very small, meaning that 
the fractional contribution of these galaxies to the total count 
of high column densities and high AVgo values is minor. If 
we assume that THINGS misses all galaxies with inclinations 
i greater than 80°, this would lead to a total incompleteness 
in the number of galaxies of cos 10° = 17% (assuming in- 
finitely thin disks). However, since the area A of a disk scales 
as A ~ cos i, the total area missed in these highly inclined 

galaxies is only cos i sin idi/ J^^ cos i sin idi ~3%. We 
conclude that the contribution from the missing highly in- 
clined galaxies to the DLA cross-section is not a significant 
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concern for this analysis. 

2.3. Beam Smearing 

If the kinematics of the H I change on spatial scales smaller 
than the beam size, measurement of the gas properties are 
subject to syste matic uncertainties. This e ffect is known as 
beam smearing. iRvan- Weber et al.l (l2005h demonstrate that 
decreasing the spatial resolution of 21 -cm data systematically 
increases the measured median AVgo. This result was estab- 
lished by convolving a high resolution H I data cube of the 
Large Magellanic Cloud (LMC) with Gaussian beams of var- 
ious widths to simulate data with different spatial resolutions. 
The increase in AVljo with beam width can be understood in 
terms of including adjacent pockets of gas that are offset in 
velocity. The resulting velocity width measurement of the to- 
tal area will be greater than the AVgo values of the individual 
Hi clouds. The LMC data has a physical resolution of 15 
pc beam^^ and a velocity resolution of 1.6 kms^^. In com- 
parison, the mean physical resolution of THINGS is 200 pc. 
Given this factor of ~ 15 difference in ph ysical resolution, 
and em ploying the same technique used by R yan-Weber et al.l 
(l2005h . we estimate that the THINGS velocity width measure- 
ments could be overstated by « 40%. Despite this expectation 
of an increase, the velocity dispersion and measurements from 
the LMC and THINGS data are similar. Specifically, the dis- 
tribution of velocity dispersion values in the LMC peaks at 
6 kms~^, wheras values derived from the THINGS galaxies 
range from 6 to 10 km s~^, i.e. not larger in all lines-of-sight. 
Thus, we conclude that beam smearing may have a mild effect 
on the THINGS data, but is unlikely to cause a major system- 
atic shift in the measured AVgo values. Of course the physical 
scales sampled by the highest resolution 2 1 -cm measurements 
are still larger than those probed by QSO lines-of-sight. In any 
case, as we will see in Section[3] if we have overestimated the 
velocity widths of galaxies at z = using the THINGS data 
this would imply an even stronger discrepancy between the 
true low redshift velocity widths and the high redshift AVgo 
values of in DLAs. 

2.4. The column density distribution function 

We have argued that THINGS provides the most suitable 
sample of 21 -cm data cubes of nearby galaxies to perform the 
kinematic analysis. However, the number of galaxies in this 
sample is an o rder of mag nitude smaller than the WHISP*" 
sample used in IZwaan et aXi (i2005bl) . which provides a very 
accurate measurement of the local H I column density distri- 
bution function /{Nm). 

To test whether the THINGS sample contains biases in the 
Hi column densities it samples, we first calculate the Hi 
column density distribution function /(A'hi)- This function 
is defined such that f {NYii)dNiiidX is the number of sys- 
tems with column density between Nm and iVni + dA^Hi 
over a distance interval dX. We use the method described by 
IZwaan et al. (2005b), where /(A'hi) is calculated as follows: 
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Here, 9(A/Hii) is the space density of galaxy i, measured 
via the Hi mass function. For the Hi mass function we 
adopt the most recent HIPASS measurement as published by 

^ The Westerbork H i survey of Spiral and Irregular galaxies 
Ivan der Hulst et ani200lh is an H I survey of 355 galaxies carried out with 
the Westerbork Synthesis Radio Telescope. 
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Fig. 1 . — The H I column density distribution function at z = from 21- 
c m emission-line obser vations. The open symbols represent the calculation 
of Zvyaa n et alj )2005b ') based on the WHISP sample, whereas the solid cir- 
cles represent the calculation presented in this paper based on THINGS data. 
Error bars include uncertainties in the H I mass function as well as counting 
statistics and indicate Icr uncertainties. 



IZwaan et"an (12005 al) . The function u'(AfHii) is a weight- 
ing function that takes into account the varying number of 
galaxies across the full range of Mm, and is calculated by 
taking the reciprocal of the number of galaxies in the range 
logA/Hii - A/2 to logMHi,; + A/2, where A is taken in 
this case to be 0.35 dex. Ai{\og A'hi) is the area function that 
describes for galaxy i the area in Mpc^ corresponding to a 
column density in the range log A'hi to log Ahi + A log Ahi. 
In practice, this area function is simply calculated by sum- 
ming for each galaxy the number of pixels in a certain log Ahi 
range multiplied by the physical area of a pixel. Finally, c/Ho 
converts the number of systems per Mpc to that per unit red- 
shift. 

Figure [T| shows the resulting column density distribution 
function /(Ahi) measured from the THING S data (sohd 
points), compared to the one published by IZwaan et alj 
(I2005bl) based on the WSRT maps of the WHISP sample 
(open points). Given the different sample definitions of 
THINGS and WHISP and the different resolutions achieved 
in both surveys, there is a remarkable agreement between the 
two data sets, reaffirming our earlier conclusion that selection 
effects are not an important issue for the THINGS sample, 
and thus that this sample is well suited for this kind of analy- 
sis. Note that we only calculated the /(Ahi) points above the 
DLA column density limit of Ahi = 2 x 10^'^ cm~^ (corre- 
sponding to log Ahi — 20.3), since this is the limit we apply 
for our kinematic analyses. In principle, the THINGS H I col- 
umn density sensitivity would allow us to calculate /(Ahi) 
down to much lower Ahi values. 

Some differences between the two samples can be seen at 
high column densities, above log Ahi ~ 21.7. This is due to 
the dearth of high inclination galaxies in THINGS. Although 
the effect shows up clearly in the logarithmic representation 
of /(Ahi), we showed in section l272l that the effect of miss- 
ing edge-on galaxies on the total H I cross section and on the 
AVgo calculation is not important for this analysis. 

The zeroth moment of /(Ahi) gives the redshift number 
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Fig. 2.— Examples of H I maps and AVgp calculations for four of the THINGS galaxies: NGC 925, NGC 3077, NGC 2841, and NGC 2903. The top row 
of panels shows the Hi column density distributions. The greyscales range from log A^jji = 19.8 to 21.8. The bottom row of panels shows the distribution of 
AVgo over the region of the galaxies for which the H I column density exceeds the DLA limit of 2 X 10^" cm~^. The full range of the greyscale is to 100 
kms^^. The black regions are those where AVgo exceeds 100 kms~^. The black spots that can be seen for example in the northern and southern outskirts 
of NGC 2903 are due to noisy spectra, where the AVgo calculation is unreliable. In NGC 2841, the extended dark regions in the AV90 map are associated with 
high Afui regions and indicate true high values of AV90. 



density dN/dz of gas above the DLA Hi column density 
threshold in the z = universe. This quantity simply rep- 
resents the number of encounters with H I gas above the DLA 
limit that a hypothetical observer would have if they were 
traveling through the local universe for a total path length 
of Az = 1. From the THINGS sample we find dN/dz = 
0.047 ± 0.009, compared to dN/dz = 0.045 ± 0.006 from 
the WHISP analysis, which again is in excellent agreement. 

2.5. The velocity widths 

The standard metric used in DLA studies to characterize the 
profile width of metal absorption lines is AVgo , the velocity 
width that encompasses the central 90% of the total optical 
depth seen in the line. This definition assures that weak ve- 
locity components that are not good tracers of the large-scale 
velocity field, as well as statistical fluctuations, are not taken 
into account. In practice, AVgo is calculated by taking the 
cumulative optical depth distribution, and finding the veloci- 
ties at which this distribution crosses the 5% and 95% levels. 
In DLA studies, the data used for this analysis are normally 
limited to those with high signal to noise, where the peak op- 
tical depth exceeds 20fT. These spectra are mostly products 
of the High Resolution Echelle Spectrometer (HIRES) on the 
Keck telescope and the Ultraviolet and Visual Echelle Spec- 
trograph (UVES) on the Very Large Telescope (VLT). This 
signal-to-noise level is calculated on smoothed spectra with 
typical velocity resolutions of 15 to 20 km s"^. 

In order to make an equitable comparison, we first smooth 
the THINGS data cubes in velocity with a boxcar filter, to a 
velocity resolution of approximately 15 to 20 kms^^, vary- 
ing slightly from cube to cube, depending on the original 
channel width. Next, we calculate normalized cumulative ve- 
locity profiles at every spatial pixel in the cube. Using the 



same strategy as defined above (just cutting the distribution at 
the 5% and 95% levels) does not give good results, because 
for the lower column densities the H I data typically possess 
lower S/N than the optical spectra (of which only the highest 
S/N examples are regarded). Therefore, the cumulative dis- 
tribution can cross these levels several times, leaving us with 
no unique definition of AVgo. To avoid this problem, we first 
identify the 50% level and travel outward in both directions 
to find the first instance where the 5% and 95% levels are 
crossed. The difference between the two velocities at which 
this happens is then defined as AVgo- In principle, this method 
could cause an underestimation of AVgo, but tests with syn- 
thetic spectra show that this effect is minimal. 

Figure ID shows four examples of the AVgo calculations for 
the THINGS galaxies. We choose to show the nearly face-on 
barred spiral NGC 925, the M81 -group member NGC 3077 
know for its extended tidal HI, the isolated nearly edge-on 
early-type spiral NGC 2841, and the barred nuclear starburst 
galaxy NGC 2903. The top row of panels shows the H I col- 
umn density distribution for these galaxies. The bottom row 
of panels shows the distribution of AVgo over the regions of 
these galaxies where the H I column densities exceed the DLA 
threshold of A'^hi = 2 x 10^°cm~^. The black regions are 
those where AVgo exceeds 100 kms^^. The black spots that 
can be seen for example in the northern and southern outskirts 
of NGC 2903 are due to noisy spectra, where the AVgo calcu- 
lation is unreliable. In NGC 2841, the extended dark regions 
in the AVgo map are associated with high A'hi regions and 
indicate true high values of AVgg. 

Examples of velocity profiles through THINGS galaxies are 
shown in Figure[3] We chose 25 random profiles, with AVgo 
values between 31 and 165 kms^^. This figure gives an idea 
of quality of the THINGS data and the accuracy of the AVgo 
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Fig. 3. — Representative Hi profiles from the THINGS data. The pro- 
files are sorted from bottom to top by increasing AVgo values. The shaded 
areas correspond to the inner 90% of the optical depth and define the mea- 
surement of AVgo- The number in the upper left corner of each panel shows 
the value of AVgo for each profile, whereas the number in parenthesis shows 
AVgo based on second crossings of the 5% and 95% levels (see text). The 
dashed vertical lines show the velocities at which the profiles reach these sec- 
ond crossings. The vertical scaling is arbitrary, the horizontal axis runs from 
-100 to 100 km s-l for all profiles. 



measurements. The shaded areas correspond to the inner 90% 
of the optical depth and define the measurement of AVgo- The 
number in the upper left corner of each panel shows the value 
of AVgo for each profile. To indicate the reliability of our 
method for determining AVgo by using the first crossing of the 
5% and 95% levels from the center, we also show in parenthe- 
sis the values of AVgo based on the second crossing of these 
levels. For most profiles, the cumulative distribution crosses 
the 5% and 95% levels only once, in which case the definition 
of AVgo is unique. In the few cases where the profiles are 
more noisy and a second crossing occurs, it is clear that the 
AVgo measurements based on the first crossings are superior; 
measurements based on the second crossing clearly overesti- 
mate AVgo- These example profiles show that our method of 
determining AVgo is accurate and not introducing a signifi- 
cant bias. 

2.6. Do low ionization metal lines trace only the neutral gas? 

Before presenting the results, we would like to come back 
to the statement made in the introduction that low ionization 
metal lines such as Fe+, Si+, Ni+ trace the neutral gas only. 
The fundamental premise of this paper is that these lines trace 
the kinematics of the high column density neutral gas, and can 
therefore be compared directly to local 21 -cm measurements. 

It is well established that approximately 25% of a spiral 
galaxy's hydrogen mass is in a layer of ioni zed gas, known 
as the diffuse ionized gas (DIG, e.g. [Reyno lds 1993). The 
scale height of this gas is 0.9 kpc in the Milky Way Galaxy, 
much larger than that of the neutral gas seen in 21 -cm. If a 
large fraction of the optical depth seen in low-ionization lines 



observed in DLAs would arise in such a medium, our results 
might be biased because this medium does not contribute to 
th e profile of the H l 21 -cm line. 

ISembach et al.l (l2000h present model calculations of ioniza- 
tion corrections in the DIG of the Milky Way, and conclude 
that elemental abundances in ionized gases might be biased 
without significant ionization corrections. However, much 
more relevant in the context of the present paper are studies 
directed at DLAs specifically since our concern is whether the 
ionized gas influences the low-ionization metal line measure- 
ments of DLA profiles. A main motivation for such studies 
of ionization effects in DLAs is the observation that Al+^ 
is present in DLAs with the same radial velocity distribu- 
tion as that of low ioni zation Unes (|Lu et al.lll996h. How- 
eve r, the results of e.g ., IViegasI (Il995h . IVladilo et all (1200 Ih 
and |Lopez_et"an (l2002l) show conclusively that ionization ef- 
fect corrections for low-ionization lines observed in DLAs are 
negligible. The direct implication of this is that the velocity 
width measured in low-ionization lines is not affected signif- 
icantly by an ionized medium. For systems below the DLA 
H I column density limit, the situa tion might be very different 
jDessauges-Zavadskv et al.ll2003h . but these are not consid- 
ered in the present analysis. 

3. RESULTS 

We first use the THINGS data to calculate the probability 
distribution function of the velocity width AVgo of galaxies 
responsible for high column density H I absorption at z = 0. 
This function defines the number of systems above a certain 
column density limit that would be encountered along a ran- 
dom 1 Mpc path through the z = universe, as a function 
of AVgo. The function is calculated by multiplying the cross 
sectional area of H I as a function of AVgo by the space den- 
sity of galaxies. If DLAs aX z = arise in galaxies similar to 
those observed in the local universe, this function would de- 
scribe the expected distribution of AVgo values for DLAs. For 
the details of calculating the pro bability distribution function 
we refer to Zwaan et al.l (l2005bl) . 

Figure |4] shows the result. The lines show the distributions 
above column density limits of logiVni = 20.3, 20.6, 20.9, 
21.2, and 21.5, from top to bottom, respectively. For all col- 
umn densities above the DLA limit of log A^hi = 20.3 the dis- 
tribution peaks at ^ 30 kms^^, and shows a long tail toward 
higher AVgo values. As one progresses to higher column den- 
sity cut-offs, the peak of the distribution shifts toward larger 
AVgo values, but the upper tail remains. We will discuss this 
further in section [TT] 

Immediately obvious from this distribution is that the peak 
in AVgo is much lower than what is typically seen in DLAs 
at higher redshifts. This is illustrated directly in Figure |5] 
where we show on a linear scale the total probability distri- 
bution function of AVgo for z = gas above the DLA limit 
(as a dashed line), together with the measured distribution of 
AVgo in DLAs at reds hifts between z = 1.5 and z — 4 from 
iProchaska et al.l (l2007 h. Note that the dashed line is the same 
as the top solid Une in Figure IH but now on a linear scale. 
The shape of the two distributions is very similar - a sharp 
peak with a high AVgo tail - but the location of the peak is 
approximately a factor two higher for the DLAs than for the 
local galaxies based on THINGS. 

This result reaffirms the conclusions of IProchaska et alJ 
(l2002h . which showed that the gas kinematics in local galax- 
ies are inconsis t ent w ith the majority of DLA sightlines. 
IProchaska et al.l (|2002|) concentrated on the gas kinematics 
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Fig. 4. — The distribution of velocity widths AV90 of high column density 
gas 2 = from the THINGS data. The lines show the number of systems as 
a function of AV90 encountered along a random 1 Mpc sightline. The right 
axis shows the corresponding number of systems per unit redshift dN/dz. 
Under the hypothesis that DLAs arise in gas disks similar to those seen in 
the local universe, these curves show the expected AV90 distributions for 
DLA absorbers. The different lines correspond to different column density 
limits: the thick, uppermost curve corresponds to the classical DLA limit of 
logTVHi > 20.3, the other curves are for log Afni > 20.6, 20.9,21.2,21.5, 
respectively. 

of the Large Magellenic Cloud and we have shown that the 
same conclusion holds if a more representative sample of lo- 
cal galaxies is considered. Whereas the kinematics of H I gas 
above the DLA limit at z = is dominated by ordered rota- 
tion in cold disks, the velocity fields of DLA systems at high 
redshifts must be influenced by other physical processes such 
as outflows, superwinds, and interactions. We discuss this in 
more detail in section |5] 

3.1. The relation between AV90 and Hi column density 

It is interesting to investigate how the velocity spread of 
the Hi gas depends on the Hi column density range con- 
sidered. The contours in Figure |6] show the probability dis- 
tribution of H I cross section in the A^hj- AV90 plane. The 
THINGS data were used to calculate the cross sectional area 
contributed by each element d^VnirfAVgo on a fine grid in 
the A^Hj-AVgo plane. The H I mass function has again been 
used to assign weights to individual galaxies. The weighting 
scheme is similar to that used for the calculation of /{Nm), 
which is expressed in Equation[T] The points in this figure are 
the DLA data from IProchaska et al.l (|2007|), over the redshift 
range z = 1.5 to z = 4. 

There are two notable effects visible in this diagram. First, 
the peak of the z ~ Q distribution occurs at values of AVqq 
around 30 km s^^, and increases slowly toward larger H I col- 
umn densities, as already seen from Figure |4] For Gaussian 
profiles, the value of AV9o=30 kms^^ translates to a veloc- 
ity dispersion of « 9 kms^^, which is a typical value for 
the gas disks in spiral galaxies. Observationally the line- 
of-sight velocity dispersion ranges from « 6 km in the 
outer parts to fti 12 km s~^ in the center (e. g., Meur er et al.l 
Il996t Ide Blok & Waltedl2006t iFraternaiieTal. 2002) . From 
the THINGS data it is found that the dispersions in the re- 
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Fig. 5. — The distribution of AV90 measured in local galaxies (dashed 
line) a nd in DLAs (histogram). The DLA data are taken from Prochaska et al] 
120071) . The local galaxy data are based on THINGS and only include sight- 
lines where the H I column density exceeds the DLA limit of A^hi = 2 x 10^*^ 
cm~^. The galaxy data are weighted such that they represent cross-section 
selected gas. The dashed line is the same as the top solid line in FigureO but 
now on a hnear scale. 



400 I 




log [cm-2] 



Fig. 6. — The bivariate distribution of AVijo and H I column density. Con- 
tours are calculated from the THINGS data, and are drawn at probability 
levels of I, 5, 1 5, 30 , 60, and 90%. The points are the DLA data from 
IProchaska et alj 120071) . Two effects are obvious from this plot: 1) the peak 
of the 2 = distribution is at AV90 ~ 30kms~^, and increases slowly 
toward larger Hi column densities; 2) large values of AVgo occur mostly in 
regions of galaxies with relatively low H I column densities. 

gions away from obvious H i shells is 11 zt 3 kms~^, aver- 
aged over the galaxy disks (iLeroy et al.ll2008h . The increase 
of AV90 toward larger A^hi values is consistent with the ob- 
servation that the Hi velocity dispersion in g as disks typi- 
cally decreases with galactoce ntric radius (e.g.. lDickev et al.l 
ll990l:lKamphuis & Sa ncisT 1993), and the fact that radial H I 
colu mn density profiles usually show a dec line toward larger 
radii (ICavatte etal.ll 19941: iBigiel et al.ll2008l) . 
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Fig. 7. — The bivariate distribution of AV90 and R, the position from the 
center of the galaxy. Contours ai'e calculated from the THINGS data, and are 
drawn at probability levels of 1, 5, 15, 30, 60, and 90%. Larger values of 
AV90 are more likely to occur at small distances from the center of a galaxy. 
Please note that the small scale structure in the 1 % contour is due to noise. 



The other notable feature in Figure|6]is that the larger values 
of AV90 are mostly associated with the lower column densi- 
ties. To understand this, we need to investigate further the 
origin of the high-end tail in the AVgo distribution. There 
are several components contributing to this tail. One of the 
most obvious contributions is that of pixels near the centers 
of galaxies, where the gradient in the rotation curve is large. 
A line of sight through these regions is likely to encounter gas 
over a large range in velocity. The steepest gradients are see n 
in earlier type spiral galaxies (e.g-. lNoordermeer et alJl20(j7l) . 
where typically the H l densities are depressed in the center 
because of the conve rsion of H l to molecular hydrogen (e.g., 
IWong & B litzll20()2h . These factors combined cause a large 
contribution of low TVhi, high AV90 points. 

To illustrate that the larger values of AV90 more likely 
to occur near the centers of galaxies we constructed Figure]?] 
which shows the probability distribution of H I cross section 
in the i?-AV9o plane, where R is the distance from the cen- 
ter of the galaxy. It is important to note that this diagram is 
also made in such a way that it represents probabilities for R- 
AV90 pairs to occur for random sight-lines through the local 
universe. Therefore, R can be regarded as an "impact param- 
eter", similar to that defined in QSO absorption line studies as 
the projected distance between the position of a background 
quasar and that of the center of a galaxy giving rise to an ab- 
sorption line. Clearly, there is an upper envelope in the distri- 
bution: the maximum AV90 seen at any value of the impact 
parameter R decreases toward larger R. 

The simple fact that the lower Nm pixels show noisier spec- 
tra causes another important contribution of high AVqq points. 
In lower signal-to-noise spectra, noise peaks more easily con- 
tribute to the central 90% flux contained in the spectra. This 
fact can be seen in Figure[3] where we showed random veloc- 
ity profiles through THINGS galaxies. For the lowest AV90 
values, the profile is clearly dominated by a single compo- 
nent. For AV90 ~ 70 kms^^, the profile typically shows two 
components, or one strong component with a lower flux tail. 



This continues to be true for even higher AV90 profiles, but 
here irregular profiles seem to get increasingly important. 

Comparing the contours and the points in Figure|6l we note 
that there is some indication for the DLA data to show the 
same trend as the THINGS local galaxies data: higher values 
of AV90 are more likely to occur at lower Hi column densi- 
ties. 

3.2. Qualitative comparison of DLA and galaxy velocity 
profiles 

So far, our analysis has concentrated on the parameter 
AV90, but it is interesting also to consider the shape of the 
profiles in more detail. In the previous section we noted that 
for sightlines where AV90 ^ 70 kms^^, the profile typi- 
cally shows two components. In fact, it appears that the in- 
dividual line components have widths within the main peak 
of the AV90 distribution in Figure 21 and the high -AV90 tail 
is mostly due to line splitting and multiple components. 

Interestingly, this same behavior is seen in the metal ab- 
sorption line profil e s of DLAs. For example. Fig 1 of 
iProchaska & Wolf3 (11997 ') shows 17 examples of velocity 
profiles of low-ionization transitions from DLAs. The nar- 
row profiles are characterized by a single component with 
AVgo ^ 30kms"^, very similar to what we observe in the 
local universe for most of the sightlines through galaxy disks 
(cf. Figure |5]l. For DLA profiles with larger widths, a shoul- 
der or tail becomes progressively more important, and a pro- 
nounced second peak starts to separate from the main compo- 
nent for AVgo values exceeding ^ 100 km s^^. Apparently, 
the second peak in Figure|5] beyond AVgo ~ 100 km s^^ is 
mostly due to profiles clearly showing multiple components, 
whereas the main peak corresponds mostly to profiles dis- 
playing a single component with a weaker tail. The profiles 
with multiple components are seen at z = as well (see 
FigureO, but are very rare. When we consider just the pro- 
files that correspond to the main peaks in Figure|5] we see 
that at 2; = they are mainly symmetric, single component 
profiles, whereas the DLAs at higher redshift typically show 
a tail. This characteristic is important for understanding the 
mechanisms repsonsible for shaping the profiles and will be 
discussed more in sectionlH 

3.3. Is AVgo good estimator of galaxy mass? 

There is accumulating evidence for the existence of a rela- 
tion between DLA (and sub-DLA) velocity spread and metal- 
licity at high r edshifts (M0ller et al. 20Q4 yLedoux et all2006l: 
iMurphy et aP 12007: Prochaska et al.. .200% . The interpreta- 
tion of this correlation, assuming that AVgo is correlated with 
mass, is that the well-known luminosity-metallicity relation 
observed in local galaxies (e.g., Tremonti et al. 2004) was al- 
ready in place at redshifts z = 4 all the way d own to z = 2. 
As a p roxy for luminosity (or galaxy mass), iLedoux et alj 
(|2006|) use measurements of AVgo. The authors argue that 
AVgo measured in low ionization lines is dominated by mo- 
tions governed by gravity, irrespective of whether the gas is 
following the rotation of the galaxy or related to infall or out- 
flow of gas or to merging of galaxy sub-clumps. In principle, 
a large collection of sightlines through the gaseous regions of 
galaxies should show a relation between depth of the potential 
well and AVgo. 

Only at z = are we in the situation that we can test 
whether the assumption that AVgo is a statistical measurement 
of galaxy mass or luminosity is valid. The result is shown in 
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Fig. 8. — The conditional probability distribution of AVgo versus absolute 
B-band magnitude, Mb- The thick solid line represents the median of the 
velocity width AVtjo distribution at each value of Mg. The thin solid lines 
show the 25th and 75th percentile distribution, the dotted lines show the 10th 
en 90th percentile, and the dashed line shows the 99th percentile distribution. 
The curves indicate global trends, small scale structure in the curves is due to 
small number statistics. The median AVgo remains constant around AVgo = 
30 kms-i. 



Figure |8] Here we show the conditional probability distri- 
bution of AVgo as a function of absolute S-band magnitude. 
This figure a created by first calculating the bivariate distribu- 
tion of cross-section in the AfejAVgo plane, and then normal- 
izing at each Mb the distribution function of cross-section as 
a function of AVgo. 

Obviously, the likelihood of finding a high AVgo value goes 
up only marginally toward brighter galaxies: the median does 
not rise significantly from low luminosity to high luminosity, 
but for more luminous galaxies the high AVgo tail becomes 
more pronounced. Therefore, in a statistical sense, AVgo is 
a very weak proxy for host galaxy mass. At any absolute 
magnitude, the median value of AVgo is ~ 30kms^^, which 
for a Gaussian profile translates to a velocity dispersion of 9 

In conclusion, at z = the velocity profile measured along 
a thin line of sight at a random position through a galaxy gas 
disk is a very weak indicator of the galaxy mass or luminosity. 
Although we have presently no means of testing this indicator 
at high redshift, the local relation does not lend much support 
to the suggestion that the observed AVgo-metallicity relation 
at high redshift can be interpreted as a mass-metallicity rela- 
tion. 

4. INTERPRETATION 

Modeling of DLA velocity profiles already suggested that 
high redshift DLAs cannot simply arise in rotationally sup- 
ported cold gas d isks such as those observed in l ocal galaxies 
(iProchaska & Wolfe 19971: iHaehnelt et al.lll998l) . "In this pa- 
per we have demonstrated empirically that the typical veloc- 
ity profiles of z = H I gas above the DLA limit have lower 
widths than those seen in DLAs. T his result is fully consistent 
with that of Prochas ka et al.l (l2002h . based on the analysis of a 
single low mass galaxy. How can these large velocity widths 
in DLAs be interpreted? 



It is well established that most actively star-forming high- 
redshift galaxies show large-scale outflows, with velocities of 
several 100 kms^^ up to 1000 kms^^. Several authors have 
suggested that these outflows could be related to a large frac- 
tion of the DLA cro s s-section a t these redshifts (N ulsen et al.l 
ll998HSchayel2001h . iNulsen et al. (1998 ) argued that the star- 
burst accompanying the formation of dwarf galaxies blows 
out most of the available gas, which results in a weakly col- 
limated wind. Once the starburst has subsided, the (largely 
ionized) ejected gas recombines and shows H I column densi- 
ties typical of DLAs. Schave (2001) built on these results and 
showed that the incidence rate of DLAs is in agreement with 
the wind scenario. He concluded that both disk gas and winds 
contribute to the high redshift DLAs cross section. 

4.1. Hi shells 

At z = 0, less extreme examples of star-formation driven 
blow-out can be found in Hi shells, which are ubiquitous 
in high resolution 21 -cm maps of nearby galaxies (see e.g. 



Brinks & Baiaia 1986; PucheeLalJ ll992riWalter & Brinksl 
X999; Kim et al. 1999). Radiative and kinematic pressure of 
stellar winds and supernovae from young massive stars in OB 
association s are thought to l ie at t he origin of these expand- 
ing shells. IProchaska et al.l ( l2005h looked specifically at the 
issue of whether AVgo measurements are affected by these 
shells in the 30 Doradus region in the LMC. It was found 
that velocity widths of sight-lines through the Hi bubbles 
were only enhanced by approximately 15 kms^^ with re- 
spect to the average LMC value. This result may not come 
as a surprise, because the typical expansion velocities of the 
shells range from 10 to 30 kms^^ in the LMC (Kimetal] 
Il999l). Similar numbers are found for other nearby galaxies 
(e.g..l Walter & Brinkslll999l) . Not only is the velocity width 
small in s hell regions, the area fil ling factor of H I shells is 
also small. 'Prochask a et alJ (l2002l) quote a covering factor of 
« 20% for regions above the DLA A'^hi limit. 

If the physics of supernova feedback at high redshifts is 
similar to what is seen in local galaxies, the above results im- 
ply that Hi shells are unlikely to explain the high velocity 
widths observed in DLAs. In order to explain the kinematics 
in high-z DLAs with outflows, one has to invoke more ex- 
treme outflow scenarios, such as that seen in star-burst galax- 
ies. 

4.2. Superwinds and tidal features 

To investigate this further, we focus on M82, the proto- 
typical starburst galaxy in the nearby universe. This galaxy 
displays a well-documented example of the superwind phe- 
nomenon, driven by the molecular driven central starburst, 
which is very prominent i n X-r ay and Ha emission (e.g., 
Martin 1998; Lehnert et aLl ll999h . M82 is also a member of 
the M8 1 group of galaxies and experiences a stron g tidal in- 
teract i on with its neighb ors M81 and NGC 3077 ( lYun et alj 
119941; IWalter et alj|2002l) . This interaction results in several 
H I tidal features in the close vicinity of M82. At a distance of 
only D — "i.Q Mpc, the galaxy can be observed in great detail, 
and therefore has been the object of studies at wavelengths 
ranging from the X-rays to the radio (see Engelbracht et alJ 
2006, for references). Here, we investigate the H I distribution 
in M82, making use of the ori ginal C - array and D-array Very 
Large Array (VLA) data of YuiTetaD (Il993h . augmented with 
B-array data from the VLA archive. The final H I column den- 
sity ni ap is very similar to the one published in iTavloret al] 
(l200l . 
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Fig. 9. — The distribution of AVgo values in M82. Only the regions above 
the DLA column density limit are shown. The central region is blanked be- 
cause of the strong radio continuum source. The full range of the greyscale is 

to 200 km ^ . The dashed line indicates the major axis of the stellar distri- 
bution in M82. Regions of AV90 « 100 kms~i and AV90 « 200 kms'^ 
are indicated in the north-east tail and in the central regions, respectively. 

We use the same procedure as described for the THINGS 
galaxies to calculate AV90 values from the M82 data for the 
regions of the H I gas above the DLA column density limit. 
Figure |9] shows the distribution of AV90 in M82. We find 
that the median AVgo from these data is « 110 kms^^. 
Very high values of AVgo ~ 200 kms^^ are seen in the 
central region over an area of several kpc^. Along the mi- 
nor axis of M82 these high velocity widths are attributed to 
interactions of M 82's gaseous medium with the superwind 
dYun et al.ll993l) . Interestingly, similarly high velocity widths 
along the minor a xis are observed in CO (Tavlor et al. 20011 

1 Walter et al.1l2002h . The tidal features in M82 also show high 
velocity widths. The south-west tail that connects up to M81 
is not well seen in our map, but the north-east tail shows AVgo 
values around 100 kms^^. 

Obviously, M82's high inclination to the line-of-sight (« 
80 deg) will contribute to the high median AVgo observed in 
this galaxy. However, the gaseous regions in this galaxy that 
are clearly identified as either tidal features or those related to 
the superwind display similarly high velocity spreads. As we 
investigated only a single starburst galaxy it is difficult to draw 
quantitative conclusions on the basis of these data. However, 
we argue that M82 clearly demonstrates that the superwind 
phenomenon and tidal interactions, which are both more fre- 
quent at higher redshifts, could go a long way in explaining 
the high AVgo values observed in DLAs. Both winds and 
tidal features could contribute to the main peak of the DLA 
distribution in Figure|5] as well as to the high AVgo bump. 

In this con t ext it is interesting to mention the results of 
iBouche et al.l (l2006l) . who find that the equivalent width of 
Mgll absorbers is anti-correlated with the mass of the halo 
that hosts them. These authors interpret this as evidence that 
the Mgll absorbers are not virialized in the gaseous haloes of 
the galaxies, but, instead, are produced by superwinds. A fair 



fraction of these Mgll systems are DLAs. iRao et al.1 (l2006h 
showed that 35% of the Mgll systems with EW (Mgll A2796, 
Fell A2600) > 0.5 A have Hi column densities above the 
DLA threshold of 2 x lO^Ocm^^. 

Additional evidence for a contribution to the DLA cross- 
section from tidal gas in galaxy groups like M81, comes from 
the recent obser vations of a much s maller sample of strong 
Mgll systems bv lNestor et al.l (|2007|) . These authors find that 
these very strong intermediate redshift Mgll absorption sys- 
tems often arise in pairs or small groups of galaxies. An es- 
timated 60% of these systems are DLAs. iNestor et all (l2007l) 
argues that the likely origin of the high equivalent width Mgll 
absorption is kinematically disturbed gas around interacting 
galaxies. 

The different peaks in the distribution of velocity width in 
the low and high redshift samples in Figure |5] could be inter- 
preted as being due to an increase in the numbe r of interact- 
ing systems with redshift. iGottlober et al.l (12001') use N-boby 
simulations to show that the number of interacting systems 
(i.e. major mergers) scales as (1 + z)^. Thus, we would ex- 
pect a significant increase in lines-of-sight with large velocity 
widths at higher redshifts. 

4.3. High velocity clouds 

iHa ehnelt et al.l (11998') have argued that irregular pro- 
togalactic clumps could (partly) produce the large ve- 
locity widths of DLAs at high redshift. Perhaps the 
most closely related examples of such clouds in the 
nearby universe are Ga lactic high velocity clouds (HVCs, 
IWakker & van Woerden|[T991i) . The H l column densities in 
the clouds themselves very rarely exceed the DLA limit, 
but it is conceivable that along a sight-line through the gas 
disk, HVCs could contribute sufficiently to the integrated op- 
tical depth to influence the AVgo measurements. Several 
authors have noted the similarity between the properties of 
Galactic HVCs an d those of Mgll ab sorption systems (e.g. 
iRichter et alll2005l;lBouche et al.ll2006l) . of which DLAs aie a 
subset. 

Because of the unknown distances to the vast majority of 
the Galactic HVCs, their physical properties are difficult to 
ascertain. In this respect, the recent identification of HVCs 
around M31 is particularly interesting (Thilker et al. 2004|). 
Synthesis folow-up observations by iWestmeier et al.. (.2005h 
show that these clouds have typical sizes of 1 kpc and column 
densities in the range 10^^ to 10^° cm~^, i.e., below the DLA 
threshold. If these HVCs really are the relics of galaxy forma- 
tion, it is likely that they present a much higher cross-section 
at higher redshifts when galaxy formation was more actively 
taking place. As such, they could have a significant effect on 
the AVgo measurements in high-z DLAs. We note that this 
component is closely related to (and perhaps partly a variant 
of) the tidal streams and winds discussed above. 

Also interesting in this context are the recent very deep 
21 -cm obse rvations of NGC 2403 (iFraternali et al.ll2002l) and 
NGC 891 (Oo sterloo elaLl l2007h . These observations re- 
veal low column density extended H I haloes, which typi- 
cally contain 10 to 30% of the total Hi mass of the galax- 
ies^. Again, the H I column densities in the haloes are too low 
to cause DLA absorption by themselves, but along sight-lines 

^ Related to these observations are the results of IHeald et al.l I200fl) who 
show evidence for diffuse ionized extraplanar gas, decoupled from the main 
rotation of the disk. However, analogs of this ionized gas in DLAs at higher 
z are unlikely to add much to the velocity widths measured in the low- 
ionization lines, for reasons explained in section l2!6l 
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through the gas disks, the haloes can have a significant ef- 
fect on AVgo measurements. These haloes show structures 
such as clouds and filaments, which resemble the proper- 
ties of Galactic HVCs. Detailed modeling of the H I haloes 
(iFraternali & Binney 2006) shows that their origin lies partly 
in a fountain mechanism ( Sha piro & Field 1976). as demon- 
strated by the close link be tween the extra-planar H l and Ha 
emission (lBoomsmall2007l and references therein). The driv- 
ing force behind this fountain is star formation in the disk. 
The cosmic star formation rate density was much higher in 
the past, up to a fact or 20 higher at z = 2 — 3 compared 
to the present epoch (iHopkins et al.ll2005h . The increase in 
the cosmic neutral hydrogen mass density is much less dra- 
matic, approximately a factor 2 over same redshifts range 
dProchaska et al.l 120051) . These two results combined imply 
that the energy input into the cold gaseous medium per unit 
gas mass was much higher in the past. This may lead to a rel- 
atively high contribution of fountain gas to the kinematics of 
cold gas at these redshifts. This high velocity gas would prob- 
ably be characterized by profiles very similar to the ones ob- 
served in DLAs displaying a main peak and a weaker tail. In 
passing, we note that the higher energy input into the gas disk 
would also lead to an increase in the gas turbulence, which 
also causes larger velocity widths. 

In order to fully explain the properties of the gaseous 
haloes, Fraternali & Binnev (2006) speculate that accretion 
from the surrounding intergalactic medium is also required. 
This "cold mode" of gas accretion is increasingly more im- 
portant at higher reds hifts in the smooth ed particle hydrody- 
namics simulations of iKeres et al.l (12005 ), indicating that this 
halo gas could play an important role in the gas kinematics of 
DLAs. 

5. CONCLUSIONS 

In this paper we have tested whether the kinematics of 
damped Lya absorbers at high redshift are consistent with the 
idea that these systems arise in the gas disks of galaxies such 
as those in the local universe. We have examined the gas kine- 
matics of local galaxies using the high quality H I 21 -cm data 
from the Hi Nearby Galaxies Survey (THINGS). To charac- 
terize the velocity widths, we adopt the parameter AVgo, the 



width that encompasses the central 90% of the total flux seen 
along any sightline through the galaxies. The same parameter 
is used for DLAs to measure the width of the low-ionization 
lines, which are believed to trace the kinematics of the neutral 
gas. 

We can use the z — THINGS data to calculate in de- 
tail how AVgo depends on H l column density, and we find 
that larger AVgo values are more likely to occur at low H I 
column densities. A similar (but weak) trend is seen for the 
DLAs, albeit at at higher values of AVgo- We also investi- 
gate how AVgo depends on impact parameter from the center 
of a galaxy and on galaxy luminosity. However, at present it 
is difficult to make a quantitative comparison between these 
calculations and the DLA data because of the very few iden- 
tifications of DLA galaxies in the literature. For a proper un- 
derstanding of the kinematics and other properties of DLAs, it 
is important that more DLA galaxies be identified at low and 
high redshifts. 

The most important result is that the median AVgo at z = 
is approximately 30 kms^^, a factor of two lower than what 
is seen in DLAs at redshifts z ^ 1.5 — 4. We discuss several 
processes that could increase the velocity widths in DLAs, 
and conclude that tidal tails, superwinds, accretion, as well as 
fountain gas are all likely to contribute to the kinematics of 
DLAs. Each of these processes are enhanced at higher red- 
shifts, where galaxy interactions are more common, super- 
winds are more frequently observed, 'cold mode' accretion is 
more efficient, and the star formation rate density is higher. 
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